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ABSTRACT 

Differential Ey-i variations ol up to ^ 0.2 mag on a scale of arcminutes across NGC 3201 are presented 
in the form of an extinction map. This map, created by calculating average Ey-i values for stars in 
small subregions of the field with respect to a fiducial region, greatly improves the appearance of the 
CMD of the cluster. We describe how we implemented this technique in detail with our data for NGC 
3201. A comparison between our map and that of the same region extracted from the COBE/DIRBE 
reddening maps published by Schlegel, Finkbeiner, & Davis (1998) (SFD) displays larger-scale similarities 
between the two maps as well as smaller-scale features which show up in our map but not in the SFD 
map. Several methods of determining an Ey-i zeropoint to add to our differential extinction map are 
presented. Isochrone fitting proved to be the most successful one, but it produces an average Ey-i for 
the cluster which is smaller than previously published values by l.Str. Finally, our results seem to 
support the statement by Arce & Goodman (1999) that the SFD maps overestimate the reddening in 
regions of high extinction. 

Subject headings: color-magnitude diagrams — dust, extinction — globular clusters: individual (NGC 
3201) — methods: data analysis — stars: fundamental parameters 



1. INTRODUCTION 

We are currently undertaking a survey of approximately 
15 Galactic globular clusters (GCs) with the aim of iden- 
tifying eclipsing binary stars around or below the turn-off- 
point by means of detecting brightness variations. These 
systems can provide masses for Population II main se- 
quence stars. High-quality photometry is essential for ob- 
taining reliable values for surface temperatures, luminosi- 
ties, ages, etc. Precise extinction determinations are a 
critical part of obtaining these results from our photome- 
try 

NGC 3201, by chance one of the first GCs we analyzed, 
is a low-latitude cluster located at 0:2000 = 10'*36™36.8'' 
and ^2000 = -46°24'40" or I = 277.2° and h = 8.6° (Har- 
ris 1996) with a retrograde orbit around the center of the 
Milky Way. Its nearby location {d — 5.2 kpc) and low 
concentration (c = 1.31) (Harris 1996) make it an attrac- 
tive target for photometric studies. Recent such studies 
include Cacciari (1984), Brewer et al. (1993), Cote et al. 
(1994), Covino & Ortolani (1997) (hereafter C&O), and 
Gonzalez & Wallerstein (1998) (hereafter GW). 

Due to its low-latitude position, the effects of differen- 
tial reddening across the face of NGC 3201 are quite sub- 
stantial. The existence of variable extinction was noted 
in practically all earlier studies of this cluster. Cacciari 
(1984), for instance, finds an irregular reddening distribu- 
tion with a variation of I^Eb-v = 0.03 mag in addition 
to a mean value of Eb-v — 0.21. GW report a range 
in Eb-v of as much as 0.1 mag across NGC 3201. Using 
spectroscopic data for 18 giants in NGC 3201, GW correct 
for the color spread by modelling the reddening in Eb-v 
as a plane in a and 5. 

Recently published dust infrared emission maps by 
Schlegel, Finkbeiner, & Davis (1998) (SFD hereafter) seem 
to indicate, however, that the dust distribution in the re- 



gion of NGC 3201 is too clumpy to be fitted by a linear 
function in position, even on the scale of arcminutes (see 
Fig. I). 

In this work, we attempt to remove the differential red- 
dening across the cluster by using our high-quality VI pho- 
tometry. Our procedure is, in principle, similar to the one 
used by Piotto et al. (1999). It aims to find the average 
AEv-i for subregions of the cluster field of approximately 
1 arcmin^ in size with respect to a fiducial region in the 
cluster where little or no differential reddening is apparent 
and where the overall Ey^i is small compared to the rest 
of the cluster. The pixel size (resolution) of this extinction 
map is approximately one fourth of the size (in area) of 
the SFD maps. Since the SFD maps are based on infrared 
emission from dust, our map will provide an independent 
check of whether, if at all, they tend to overestimate red- 
dening, as was suggested by, e.g., Arce & Goodman (1999) 
for regions with Ay > 0.5 mag. 

Our results indicate variations in Ey-i of up to ~ 
0.2 mag. Differential reddening this strong can wreak 
havoc with photometric and spectroscopic studies of clus- 
ter stars. The inherent uncertainties in the various param- 
eters of the CMD (e.g., magnitude and color of the main- 
sequence turnoff (MSTO), horizontal branch features, etc.) 
are greatly amplified. Moreover, one may "detect" age or 
surface-temperature gradients where there are none in re- 
ality. To give an example of the magnitude of the effect, 
using the recently published color-temperature relations 
by Houdashelt, Bell, & Sweigart (2000), the effective tem- 
perature of a solar-metallicity, main sequence star with 
V — I ^ 0.8 would vary by approximately 600K for a dif- 
ferential reddening effect of Ey-i ~ 0.2 mag, which, in 
turn, could lead to errors in the metallicity determination. 

Using our internal dereddening technique, we obtain a 
high-quality, deep-photometry CMD of NGC 3201 com- 
prised of approximately ninety 600s exposures as well as 
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some shorter ones, all using the V and / bands. 

Our observations and data reductions, as well as the 
details about our internal dereddening method, are docu- 
mented in Section 2. Section 3 contains our results con- 
cerning the reddening map and the CMD of NGC 3201. 
Finally, we discuss the determination of the reddening ze- 
ropoint of our extinction map in Section 4 and give a brief 
summary of our work in Section 5. 



2. OBSERVATIONS, 



DATA REDUCTION, 
DEREDDENING 



AND INTERNAL 



2.1. Observations and Basic Data Reductions 

The NGC 3201 observations were obtained during the 
nights of April 26 through May 6, 1998, at the Las- 
Campanas-Observatory (LCO) Im Swope Telescope, using 
Johnson-Cousins VI filters and a SITe 1 2048^ CCD with 
a field-of-view of 23.5 arcmin on a side. Table 1 gives the 
number of epochs we observed for different exposure times 
which were chosen in order to cover a larger magnitude 
range in the CMD. 

The initial processing of the raw CCD images was done 
with the routines in the IRAF^ CCDPROC package. For 
each night, 10 bias frames were combined for the bias sub- 
traction. The V band flats were produced by combin- 
ing between 4 and 6 twilight flat images per night. All / 
band data were first flattened using an image comprised 
of 15 individual / band domeflats, and then divided by a 
normalized, dome-flattened dark-sky-flat which itself was 
created by median-averaging approximately 40 individual 
different blank sky fields. Finally, the / band images were 
corrected for fringing by subtracting a fringe image which 
was created by subtracting the mean pixel value of the 
dark-sky-flat from the dark-sky-flat image itself. 

The processed data were reduced using DoPHOT 
(Schechter, Mateo, and Saha 1993) in "fixed-position 
mode" where the positions of stars are fixed (after cor- 
recting for global frame-to- frame shifts as well as small dis- 
tortions) to the positions measured on a deep-photometry 
template image obtained by co-adding the ~ 15 best-seeing 
frames. 

The photometric results for every star were averaged 
over all frames with the same exposure time to obtain the 
final magnitude for the star under investigation. For the 
600s exposures, our only requirement was that a star ap- 
pear in more than 75% of the epochs. The 60s images 
were taken to complete the CMD in the brighter regions; 
we only took magnitudes from these exposures of stars 
which were saturated in the 600s exposures. The same 
procedure was followed for the 10s frames. Aperture cor- 
rections were applied by calculating one constant value for 
the whole chip after we found little dependence upon po- 
sition on the CCD. The data from the shorter exposures 
were shifted to the photometric system of the 600s expo- 2. 
sures using non-saturated stars in common. 

During our consecutive two photometric nights (May 
2 & 3, 1998), we obtained a total of more than 140 ob- 
servations of various Landolt (1992) standard stars which 
covered a range of —0.2 < V — I < 2.2 in color as well as 
1.15 < X < 1.55 in airmass. Using the IRAF PHOTCAL 
package, we applied one single standard star solution for 

^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the NSF. 



data from both nights of the form 

V = V + Of) + boXy + co{v - i) (1) 

V - I ^ai+biX^ + ci{v -i), (2) 

where the Cj, 6^, and Cj are the fitted constants, Xfuter 
is the airmass of the exposure taken with the respective 
filter, the lowercase magnitudes are instrumental, and the 
uppercase ones are the known magnitudes of the standard 
stars. The root mean square error was 0.007 mag and 
0.016 mag for V and V — I, respectively 

Astrometry was performed by identifying 59 USNO ref- 
erence stars (Monet et al. 1996) in NGC 3201 and using 
the IRAF IMAGES package for the coordinate transfor- 
mation. A linear fit in x and y produced errors around 
0.2 arcsec, consistent with the USNO precision. The tan- 
gent point of the transformation was (x, y) — (962.7, 934.2) 
and corresponds to a2ooo = 10''17'"31.2'* and ^2000 = 
46°24 11 . The rotation of x and y with respect to a 
and 5 was 91.23°, and the pixel scale came out to be 0.695 
pix/ arcsec. These results enabled us to create the refer- 
ence grid in Fig. ||. 

2.2. Mapping the Differential Reddening 

At this point in the data reduction we were in the pos- 
session of a fully calibrated CMD whose main sequence 
appeared exceptionally broad (see Fig. ^). After ruling 
out a number of possible instrumental causes for this ef- 
fect, and after we noted a clear dependence of the appear- 
ance of the main sequence upon the positions of stars on 
the CCD (illustrated in Fig. |l|), we concluded that the 
underlying cause for the broadness of our main sequence 
was differential extinction across the field of view. This 
seemed especially sensible given the low-latitude location 
of NGC 3201 and motivated us to create a reddening map 
since our binary star data were not useful without this 
correction. In order to create such a map for NGC 3201, 
we developed the steps described below (see also Fig. H): 

1. A fiducial region was chosen in which 

• very little or no differential reddening occurred, 
i.e., where the main sequence appeared narrow; 

• the overall Ey-i was very small with respect 
to the rest of the field of view; 

• there were enough stars in the field to obtain a 
statistically reliable fit to their positions on the 
CMD. 

The coordinates of this fiducial region are approx- 
imately 10'' 16^26^ < a20oo < lO'' 16^53^ and 
-46°3l'2l" < (52000 < -46°26'44" (see Figures | 
through 0). 



A high-order polynomial was fit to the main sequence 
stars with 18 < F < 21 and 0.7 < F - / < 1.4 in the 
fiducial region, using algorithms supplied by Robbie 
Dohm-Palmer (2000, private communication) which 
enabled us to separate the main sequence stars from 
background stars in the CMD, and by Press et al. 
(1992), which provided the least-squares fit to the 
data (see Fig. 
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We divided the field of view into subregions of differ- 
ent sizes in order to have a sufiicient number of stars 
in each region. Towards the outer parts of the field, 
the stellar density decreases. As a result, the sizes 
of the subregions increase. The inner subregions are 
all of the same size (100x100 pixel). See Fig. ||. 

For each of these regions, every star falling between 
17.9 < V < 21.1 and 0.65 < V - I < 1.45 was 
incrementally moved along the reddening vector de- 
fined by the relations in Cardelli, Clayton, & Mathis 
(1989) until it intersected the fit (see Fig. |). 

The statistical biweight (see Beers, Flynn, & Geb- 
hardt 1990) of all these incremental shifts of the 
stars in a given subregion was then calculated, out- 
liers (0.9 a) were removed, and it was recalculated. 
This provided the value of the total shift along the 
reddening vector for each star. Since the slope of 
the reddening vector for a standard extinction law 
is known (1.919; see Cardelli, Clayton, & Mathis 
(1989)), the differential Ey-i for each star corre- 
sponds to the {V — /)-component of the vector de- 
scribed above (Fig. ^). 

For six of the 100x100 pixel subregions as well as 
one 200x300 pixel subregion on the extinction map 
shown in Fig. ^ too few stars were present to cal- 
culate the average reddening. For these cases, we 
calculated Ey-i by simply averaging the reddening 
values of all the neighboring subregions^. 

Our error analysis for the shifts described above in- 
volved the following: 

• We created an error ellipse for every star de- 
fined by the values of the associated errors (as 
returned by DoPHOT) in color and magnitude. 

• Since the color and magnitude errors are corre- 
lated, the error ellipse is tilted, with the tilt 
angle and the lengths of the semimajor and 
semiminor axes functions of ay a-nd ay-i- 

• One average tilt angle was obtained for every 
subregion under investigation by averaging the 
results for the individual stars. For every sub- 
region, this angle was approximately 35° ± 4°. 
The semimajor and semiminor axes for every 
subregion were obtained in the same way, so 
that every star in a given subregion has the 
same error ellipse associated with it in the cen- 
ter. 

• This error ellipse was shifted along with its re- 
spective star during the dereddening process de- 
scribed above in items 4 and 5. The "point 
of first contact" (pfc), i.e, the point at which 
the error ellipse first touches the fit through the 
data in the fiducial region, as well as the "point 
of last contact" (pic), represent the respective 
one-sigma-deviation points (see Fig. ||) . 



• These two contact points are not necessarily 
symmetric about the reddening value of the 
star (center of error ellipse), but in order to 
increase the readability of this publication, we 
calculated the mean pfc and pic for a given sub- 
region, and then their average distance from 
the center of the ellipse, so that Ueiupse = 
0.5{pfc — pic). 

• The final error estimate for each subregion 
(cfEv-i) given as the lower number in each of 
the pixels in Fig. ^ is then obtained by adding 
in quadrature aeinpse and the error in the mean 
of the shifts of all the stars in the subregion un- 
der investigation. 

3. RESULTS 

3.1. Extinction Map and Comparison with SFD Map 

Fig. H summarizes the results of our internal deredden- 
ing procedure. Darker regions correspond to higher red- 
dening with respect to the average reddening taking place 
in the fiducial region described in Section 2.2, item 1. For 
reference purposes, we included the approximate locations 
of this fiducial region and the core radius of NGC 3201, as 
well as a coordinate grid. The average Ey-i relative to 
the fiducial region is 86 mmag with a standard deviation 
of 61 mmag. The individual subregions' Ey-i values in 
millimagnitudes (relative to the average reddening in the 
fiducial region) are shown as the top number in each of the 
pixels in the grid of Fig. ^. The bottom number in each 
pixel corresponds to the error in the mean for the corre- 
sponding Ev-i value. In addition, we included an image 
of NGC 3201 (Fig. ^ in order to indicate the locations 
of the fiducial region and the core radius of the globular 
cluster. 

At this point, these reddening values are all relative to 
our fiducial region which itself is suffering some mean in- 
terstellar extinction. The procedure described cannot de- 
termine the absolute reddening zeropoint. Thus, we had 
to introduce additional information, derived from previous 
studies, in order to calculate this zeropoint. The discus- 
sion of this process is given in Section 4.1. 

As a first step, we subtracted our extinction map from 
the SFD map of the region of NGC 3201 which we ob- 
served. The SFD map of the cluster region is shown in 
Fig. ^ with the same orientation and size as Fig. ^. Once 
again, dark corresponds to higher extinction. We added 
the location of the fiducial region for reference. The aver- 
age Ev-i of the SFD map is 487 mmag with a standard 
deviation of 59 mmag. 

A first comparison indicates that, while the larger-scale 
structure is the same for both extinction maps (in partic- 
ular the obscuration ridge extending from the top left part 
of the image to the right center part of it), our map shows 
smaller-scale features which are not present in the SFD 
maps. These features are visible in the difference image in 
Fig. whose orientation and size is once again identical 
to the other Figures. Darker regions correspond to places 
where our data indicate the presence of dust which does 
not show up in the SFD map; lighter regions show good 



In this process, every 100 x 100 subregion was treated individually, even if it was a part of a larger one. That is, for the 200x300 region just 
south of the north-east corner of the field of view in which the number of stars was too low to obtain a reddening value, the above procedure 
was applied six times. 
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agreement between the two maps. The average Ey-i of 
this difference map is 401 mmag with a standard devia- 
tion of 46 mmag (significantly smaller than the standard 
deviation of both our extinction map and the SFD map). 
This lower standard deviation is an indication that both 
our map and the SFD map of the region of NGC 3201 
are tracing many of the same features in the foreground 
extinction along the line of sight to the globular cluster. 

The smaller-scale features present in only our map but 
absent in the SFD map are certainly worth further atten- 
tion. At this point, we are not entirely confident that we 
can give a definite reason for the discrepancies, but some 
potential reasons could be the following: 

• The smaller-scale features which are not in the SFD 
maps might have been smoothed out by the lower 
resolution of the IR observations used to create the 
maps. 

• The temperature of the dust producing the smaller- 
scale features was outside the detection range of the 
instrument used to collect the data for the SFD 
maps. This way, they would still show up in our 
analysis. 

• The dust we are seeing is actually part of the globu- 
lar cluster itself"^ and was therefore not picked up in 
the SFD maps. 

3.2. Photometry and CMD 

Fig. contains the raw CMD (no dereddening applied) 
for NGC 3201. To illustrate the effects of our deredden- 
ing method, we show, in Fig. ||, the CMD after apply- 
ing the extinction map (differential reddening correction) 
to the data. The datapoints are now all shifted to the 
CMD-location of the stars in the fiducial region, that is, 
no absolute reddening zeropoint has been applied. The 
improvement is quite stunning. The width of the main se- 
quence of Fig. ^has, by applying the differential reddening 
map, decreased to a fraction of its former value in Fig. pi 
The subgiant branch which was hardly visible at all in the 
raw CMD is now well traced out by the datapoints. Even 
the horizontal branch, though sparsely populated in both 
CMDs due to saturation effects of the bright stars in the 
field, is narrower in the differentially dereddened CMD. 

4. THE REDDENING ZEROPOINT OF OUR EXTINCTION 
MAP FOR NGC 3201 

4.1. Comparing the MSTO Color with Previous Studies 

In order to make our reddening map (Fig. |^) a useful 
tool, we need to determine the reddening zeropoint to add 
to the Ev~i values in the grid. As mentioned above, this 
could not be done using our data alone. In order to be able 
to use results from previous studies for the calculation of 
this zeropoint, an agreement between our photometry re- 
sults and the corresponding literature values is necessary. 
When comparing MSTO color, however, one needs to con- 
sider that different parts of the field of view are reddened 
by different amounts, as we showed in our extinction map 
(Fig. I). 

We estimated our main-sequence turnoff Vto magnitude 
to be 18.2 and the MSTO color {V - I)to to be around 

^The phenomenon of dust inside globular clusters was discussed by, 



0.88 mag; both were determined by eye from the raw pho- 
tometry CMD (Fig. I). 

These values were compared to the results of Rosen- 
berg et al. (2000) (hereafter RB) and C&O. RB find 
Vto ~ 18.2 and {V - I)to ~ 0.92; both values are 
again determined by eye. C&O report Vto ~ 18.15 and 
{B — V)to ^ 0.65 (tabulated). Using their reddening 
estimate of Eb-v ~ 0.22, the relations in Cardelli, Clay- 
ton, & Mathis (1989), and a set of isochroncs provided by 
Don VandenBerg for [Fe/H] = -1.41 (VandenBerg 2000, 
private communication, based on evolutionary models by 
VandenBerg et al. 2000; hereafter VDB), the C&O value 
for {V — I)to can be calculated to be approximately 0.93. 

One discrepancy between our results and these litera- 
ture values is therefore that our value for (V — I)to is 
approximately 0.05 mag bluer than what is quoted by the 
authors mentioned above. The existence of variable red- 
dening may be a likely explanation for this discrepancy, 
especially when one considers which part of the cluster 
was observed. 

RB concentrate on the southern and eastern parts of the 
cluster for their longer exposures and the center of NGC 
3201 for their shorter exposures. All of these three regions 
are, according to our reddening map (Fig. ^), regions of 
higher extinction with respect to our fiducial region. 

C&O observe regions just to the north and west of the 
cluster center and the cluster center itself. While the north 
and west regions suffer only low reddening (Fig. ||), the 
area within ~ 7 arcmin to the north and west of the clus- 
ter center is obscured by a smaller feature not visible in 
the SFD map (see Fig. |). At least 5 of 7 C&O fields are 
very likely affected by this dense obscuration feature. 

While a quantitative (star-by-star) magnitude compar- 
ison between our data and the data obtained by RB and 
C&O is not possible, we conclude that their estimates for 
the MSTO color could very well be influenced by the fact 
that the regions they observed were subject to extinction 
somewhat higher than the average value over our field of 
view. 

We note that these comparisons focus on the MSTO 
colors; the vertical nature of the MSTO makes it difficult 
to compare the V mag of the MSTO at sufficiently high 
precision. 

4.2. Basic Isochrone Fitting and Determination of 
Reddening Zeropoint 

We examined three different possibilities of finding the 
Ev-i zeropoint to add to our extinction map. All three 
methods arc outlined below. For each case, it is worth 
pointing out that a direct comparison with literature val- 
ues should be taken with caution since usually, only a sin- 
gle numerical value for the reddening is given, whereas our 
result is a reddening matrix of which every single element 
needs to be added to a zeropoint offset. 

4.2.1. Using Isochrone Fits 

Simultaneously fitting age, distance, and reddening to 
the aforementioned set of VDB isochrones produced the fit 
shown in Fig. |l^. It should be noted that since the point 
of this publication is outlining our dereddening method as 

e.g., Forte et al. (1992). 
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well as creating an extinction map for NGC 3201, we only 
performed some basic isochrone fitting here. For an esti- 
mate of [Fe/H]^ we adopted the isochrone with metallicity 
closest to the GW value of -1.42. 

For the fit of Fig. |l^, [Fe/H] = -1.41, d ~ 4.5 kpc, 
age =18 Gyrs, and the Ey-i zeropoint to be added to 
the values in Fig. | is 0.15 mag. The age is clearly on the 
high end of things, but the 18 Gyr isochrones produced by 
far the best fits, even when [Fe/H] and the distance were 
slightly modified. The value for the distance is slightly 
below the one quoted by Harris (1996), but fairly close to 
the two values derived using slightly different methods by 
C&O. 

When adding the mean Ey_i of our reddening map (86 
mmag) to the zeropoint, the average Ey-i for the cluster 
is approximately 0.24, which is slightly below the estimates 
of, e.g., Cacciari (1984) whose Ey-i ~ 0.32, and of Harris 
(1996) and RB who quote the same value. 

In order to justify our choice of isochrone parameters 
and to eliminate the possibility of multiple factors con- 
spiring against us in the attempt to find a good isochrone 
fit to the data, we illustrate in Fig. |ll| what the CMD 
would look like if we were to strictly adopt literature val- 
ues from Harris (1996) {[Fe/H] = -1.58, Ey-i = 0.325, 
d = 5.2 kpc). Plotted are a pair of isochrones with re- 
spective ages 12 and 14 Gyrs on top of the data^. The 
agreement between data and fit is not as good as in Fig. 
10, 



4.2.2. Using SFD Maps 

As stated in Section 3.1, the mean Ey-i of the difference 
map between the SFD and our data is 401 ±46 mmag. The 
average Ey-i value for this difference map corresponds to 
the Ev-i zeropoint to add to our data to obtain absolute 
reddening values. 

Calculating the total reddening in this fashion gives 
Ev-i ^ 0.49 which is significantly higher than the litera- 
ture values. Furthermore, using the isochrones mentioned 
in the previous subsection, we were not able to produce a 
fit to the data due to the fact that Ey-i was too high. 

4.2.3. Using RRLyrae Stars 

Another method of attempting to find a reddening zero- 
point is to simply adopt one from the literature and model 
the thus dereddened data with a set of isochrones. Here, 
we took Cacciari's (1984) value of Ey-i ^ 0.32 which she 
calculated using RRLyrae stars. Since this corresponds to 
the average reddening for the whole cluster, we had to sub- 
tract the mean Ey-i of her RRLyrae stars (in our field 
of view) with respect to our fiducial region. This value 
turned out to be 79 mmag. 

We failed, however, to find a good isochrone fit to the 
thus dereddened data, as is shown in Fig. |lj. In order 
to get the MSTO points of the data and the isochrone to 
agree, we had to lower [Fe/H] to -2.01, a value quite a bit 
below literature values such as -1.41 (GW) or -1.58 (Harris 
1996). The distance modulus here is 13.4, and the age 14 
Gyrs. As one can easily see, the fit fails to hit the subgiant 
branch due to the "young" isochrone age of 14 Gyrs. Older 
age isochrones could not be shifted far enough toward the 

^In order to have a valid basis for comparison, we needed to apply 
data, where 0.325 is the Harris (1996) Ey_j value and 0.086 is the 



dereddened MSTO location to produce a fit, even when 
the metallicity was lowered to the lowest available value of 
-2.31. 

4.2.4. Comments on the Adopted Reddening Zeropoint 

Given that we were not able to produce a fit to the data 
using the VDB isochrones when we calculated the redden- 
ing zeropoint using either the SFD maps or the Cacciari 
(1984) RRLyrae stars, we conclude that our best estimate 
for the Ey-i zeropoint is 0.15 mag, giving an average 
Ev-i ~ 0.24 for the cluster as a whole. This value falls 
approximately 1.5(7 below the estimate by Cacciari (1984) 
and the value tabulated by Harris (1996). 

We could force our reddening zeropoint to agree with 
these past results by applying a systematic shift to our / 
magnitudes of 0.08 to brighter values which would give a 
distance to NGC 3201 of 4.8 kpc. We have, however, no 
reason to adopt such a shift. The transformations to stan- 
dard magnitudes as well as the application of the aperture 
corrections (see Section 2.1) underwent repeated thorough 
examinations and appear to be correct. Since no other in- 
dication throughout our data reduction and analysis sug- 
gests a systematic error in the / band magnitudes, either, 
we simply cannot justify adjusting our / magnitudes in 
order to compensate for the 1.5cr effect described above. 

5. SUMMARY AND CONCLUDING REMARKS 

In the process of finding eclipsing binary stars in NGC 
3201, we noticed the existence of variable reddening of up 
to 0.2 mag in Ey-i on a scale of arcminutes. Using our 
internal dereddening method outlined in Section 2.2, we 
created an extinction map which is shown in Figures I and 
m Applying the map to our raw data (Fig. ||) significantly 
improved the appearance of the CMD (Fig. H). 

Comparison between our extinction map and the SFD 
map of the same region (Fig. ^) showed that the same 
larger-scale features exist in both maps. Our map displays 
some additional, smaller-scale features which are absent in 
the SFD maps (see Fig. 0). 

The Ev-i zeropoint which needs to be added to the 
numbers in Fig. ^ in order to get absolute Ey-i values 
is 0.15. This value is below literature results (Cacciari 
1984, Harris 1996) by approximately 1.5tT, but produced 
by far the best VDB isochrone fit to the data. The ze- 
ropoint determined with the help of the SFD maps gives 
Ev-i ~ 0.49 as the average value across NGC 3201 which 
is higher than literature values which supports the state- 
ment by Arce & Goodman (1999) that the SFD maps over- 
estimate the reddening in regions of high extinction. 

The results from this work will be essential in our binary 
star research where high-quality photometry of every bi- 
nary system is necessary for distance determinations. A vi- 
tal condition to obtaining these measurements is of course 
the knowledge of the exact extinction the star under in- 
vestigation is suffering. 

Furthermore, studies like this may be useful in deter- 
mining properties of the ISM itself such as examining a 
possible dependence of Ry upon position in the field of 
view or giving insight into the distribution and properties 
of the dust along the line of sight. 

an Ey—i zeropoint of 0.325 — 0.086 = 0.239 to the differentially dereddened 
mean Ey_j value of our extinction map. 
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Fig. 1. — "Raw" (calibrated but before dereddening) CMDs of 9 subregions of the field of view of NGC 3201. The x and y 
ranges of the respective subregion are given on the axes of each individual CMD. a: = 0, 700, 1400, 2100 approximately correspond 
to (52000 = -46°13'03",-46°2l'09",-46°29'l6", and -46°37'22", respectively; y = 0,700,1400,2100 approximately correspond to 
"2000 = 10'n6™27^10'n7™14'',10''18'"01^ and 10'"18'"48^ respectively. The directional arrows on the bottom and the left of the im- 
age show the approximate orientation of the CCD. The varying appearance and broadness of the main sequence as a function of position 
clearly indicates the differential reddening across the field of view. The line through the CMDs represents the best fit through the datapoints 
in the fiducial region 1200 < x < 1600 and < j; < 400 (described in Section 2.2, item 1). It is apparent that there is not only less differential 
reddening in that region (tight main sequence), but also that the overall Ey_j is low. The reddening vector for Ey_i = 0.2 is shown in the 
top right panel. 
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Fig. 2. — This figure illustrates the dereddening procedure which was used in this paper. The triangles represent stars in the subregion 

700 < X < 800 and 800 < y < 900 (approximately 10''17'"21'' < 02000 < 10''17'"28'', as well as -46°2l'07" < ^2000 < -46°22'l6") before 
being dereddened; the filled circles represent the same stars after the dereddening along the reddening vector shown in the top right corner of 
the figure. The x-shaped symbols are the stars in the fiducial region (mentioned above) in which none or very little differential reddening is 
taking place. The fitted polynomial with range 18 < mag < 21 is visible as the curve through the fiducial region datapoints. The reddening 
for this region is close to Ey—i = 0.2 with respect to the average of the fiducial region. The tilted error ellipses, described in Section 2.2, 
item 7, greatly exaggerated in size, are shown to illustrate our dereddening method which is schematically outlined in the box to the upper 
right. The straight line represents the fit through the datapoints in the fiducial region. Example a represents the "point of last contact" 
of the error ellipse with the fit, example b the shift from the original datapoint until it intersects the fit, and example c the "point of first 
contact" between the error ellipse and the fit. 
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LO''13^]OJOP 



_ LO''13^00JOP 
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□« f3000) 
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Fig. 3. — Jhe extinction map. East is up, north is to the left. The darker the color of a subregion, the more extinction is occurring in 
it (see Fig. M). Note the "ring" of larger subregions (with lower density of stars) around the center of the field of view. The size of every 
small subregion is approximately 70 arcsec on a side. For reference, we included an illustration of the core radius of NGC 3201 (from Harris 
1996) around the location of the cluster center on the extinction map as well as the location of the fiducial region describedJn Section 2.2. 
Furthermore, we overlaid a grid for coordinate reference. The coordinate axes x and y correspond to the ones used in Fig. kj. The average 
reddening with respect to the extinction occurring in the fiducial region is Ey-i = 86 it 61 mmag. 
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Fig. 4. — The Ey—i values with respect to the fiducial region for the individual pixels of our reddening map in millimagnitudes. The pixels 
correspond to the subregions in Fig. H with the respective values for x and y coordinates. The top number in every pixel is the value for 
Ev—i obtained with the method outlined in Section 2.2. The bottom number is the associated error, as described in Section 2.2, item 7. The 
pixels where the error value is N/A are the ones with an insufficient number jof^tars, so the reddening results were obtained by interpolation 
from neighboring values (see Section 2.2, item 6). For reference to figures H, pi H, and vA (same orientation), the location of the fiducial region 
is indicated in the lower right part of the figure. The average reddening with respect to the mean reddening occurring in the fiducial region is 
Ey—i = 86 ± 61 mmag. In order to obtain a map of absolute reddening values, an Ev—i zeropoint has to be added to the data. We discuss 
the various zeropoints we examined in Section 4.1. Finally, to obtain an absolute value for Eg_v, one needs to multiply the absolute (i.e., 
including the zeropoint) Ey—j value by 0.6468 (Cardelli, Clayton, & Mathis 1989). 
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Fig. 6. — A graphical representation of the SFD data in the region ofJ\F 
given for reference; the orientation of the map is the same as in Figures H, 
ridge of obscuring material extending from the top left part of the map to 
not show up on the SFD map (see Fig. VV). The average reddening is Ev—i 



that we observed. The location of our fiducial region is 
and [? . Common features are easily recognizable such as the 
ilie'right center. A number of smaller-scale features, however, do 
= 487 ± 59 mmag. 
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Fig. 7. — A difference plot between our extinction map (Fig. bh and the SFD map (Fig. rj) with the same orientation. Darker regions 
correspond to areas where our map indicates more obscuration relative to the average value than the SFD data; lighter regions are places 
where the two maps agree very well. The average reddening is Ey—i = 401 ± 46 mmag, which corresponds to the SFD-map-zeropoint (see 
Section 4.2.2). Note that the standard deviation is lower than the one of the SFD map and the one of our extinction map. 
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Fig. 9. — The internally dereddened CMD of NGC 3201. V and V — I indicate the locatixm of the datapoints after the differential reddening 
with respect to the fiducial region was corrected for, i.e., after the inverse values of Fig. H were applied to the data. No reddening zeropoint 
is applied to the data in this plot. 
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Fig. 12. — The fit to the data dereddened by adopting the Ev-i zeropoint of 0.24, using Cacciari's (1984) RRLyraes, and applying our 
reddening map. [Fe/H^ = —2.01 (much lower than literature values), age = 14 Gyrs, and distance modulus = 13.4. As one can see, the fit 
fails to reproduce the loci of the datapoints. We therefore conclude that the reddening zeropoint of Fig. hCj is correct. 
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